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ABSTRACT 

The nature of the cosmic dark matter is unknown. The most compelling hypothesis is that dark matter 
consists of weakly interacting massive particles (WIMPs) in the 100 GeV mass range. Such particles would 
annihilate in the galactic halo, producing high-energy gamma rays which might be detectable in gamma ray 
telescopes such as the GLAST satellite. We investigate the ability of GLAST to distinguish between WIMP 
annihilation sources and astrophysical sources. Focusing on the galactic satellite halos predicted by the cold 
dark matter model, we find that the WIMP gamma-ray spectrum is nearly unique; separation of the brightest 
WIMP sources from known source classes can be done in a convincing way by including spectral and spatial 
information. Candidate WIMP sources can be further studied with Imaging Atmospheric Cerenkov Telescopes. 
Finally, Large Hadron Collider data might have a crucial impact on the study of galactic dark matter. 
Subject headings: dark matter — elementary particles — Galaxy: halo — gamma rays: theory 



1. INTRODUCTION 

It is now firmly established that the majority of matter 
in the universe is non-baryonic. Evidence for this standard 
cosmology includes the microwave background anisotropics 
(Spergel et al. 2006) and th e power spectrum of density fluc- 
tuations on galactic scales dTegmark et al.ll2006l) . The "dark 
matter" is of unknown composition, but indirect evidence 
from particle physics and cosmology indicates that it is likely 
to consist of weakly interacting massive particles (WIMPs) 
in the mass range 30 GeV to 3 TeV. Such particles would 
be expected to annihilate slowly in galactic halos. In most 
WIMP models, a large fraction of the annihilation radiation 
is expected to be gamma rays from the decays of the ir° me- 
son, produced copiously in any energetic interaction involving 
hadrons. 

In the cold dark m atter (CDM) paradigm (Blu menthal et all 
Il984t [Peebles] Q 984), it is well known that structure forms 
hierarchically: the dark halos of galaxies such as the Milky 
Way are expected to contain large numbers of sub-halos. For 
WIMPs, the sub-h alo mass spectrum is expected to exten d 
down to 10" 6 M Q dGreen et all [20051: iDiemand et alj 120051) . 
The substructure is expected to be nearly isotropic, thus anni- 
hilation in the sub-halos can be away from the galactic plane, 
where astrophysical sources are concentrated. 

The brightest source of WIMP annihilation radiation is ex- 
pected to be the galactic center (GC), where the WIMPs are 
most concentrated. The HESS collaboration has concluded 
that dark matter annihilation radiation is at most a small frac- 
tion of the emission c oming from the GC above 100 GeV 
dAharonian et alJl2006T). The MAGIC collaboration has con- 
firmed these results dAlbert et ai1l2006l) . To explain these data 
in terms of dark matter annihilation, WIMPs of order 10 TeV 
mass would be required. We will proceed with the plausible 
assumption that the GC emission is predominantly astrophys- 
ical. In fact it has been shown that there is only a narrow 
window in which dark matter annihilation could be discov- 
ered at the GC, given the HESS source and its extrapolation 
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below 100 GeV (Zahariias & Hooper 2006). This emphasizes 
the point that galactic satellites may be the most promising 
sources of WIMP annihilation radiation. 3 

Many authors have discussed the pos sibility of an- 
nihila tio ns in ga lac tic s ubs tructure ( Bergstro m et alJ 
[1999|: iBaltz et all l2000t ICalcaneo-Roldan & Moorel 
bood iTasitsiomi & Olintol 120021: IStoehr et all 1 2003 1 : 
Taylor & Silkl 120031: lEvans et all [20041: lAloisio et"aT] 12004 
iKoushiappas et alJl2004t lBil l2006t IDiemand et all 120071) : In 
this Letter, we will illustrate that WIMP annihilation sources 
are distinguishable from all known (observed) astrophysical 
source classes. The detection of a steady, extended high 
latitude source with a WIMP annihilation spectrum (e.g. by 
the GLAST satellite) would provide strong evidence that the 
dark matter in the Galaxy actually consists of particles in 
the 100 GeV mass range, a crucial piece of the dark matter 
puzzle. 

2. gamma rays from hadronic interactions 

Collisions of relativistic hadrons, such as when cosmic ray 
protons impinge on the interstellar medium (ISM), are typi- 
cally inelastic. Energy is lost mostly to 7r + ,7r°,7r" mesons, in 
roughly equal numbers. The decays of the ir° mesons domi- 
nate the gamma rays from hadronic interactions. The ir° has a 
mass m-n = 135.0 MeV, and it has two common decay modes: 
7T° -> 27 (98.8%) and ir° -> e + e~7 (1.2%) with rare modes 
contributing less than 0.01%. As a pseudoscalar particle, it 
decays isotropically. The photons emitted in 7r° — * 27 have 
energies of £0 = '«tt/2 = 67.5 MeV in the ir° rest frame. In 
the lab frame, the energies are E± = £07 (1 ± (3 cos 8cm)- The 
isotropy of the decay implies that cos 9cm is uniformly dis- 
tributed, and thus the spectrum dN/dE is constant between 
the minimum and maximum energies £ m i n max = £07(1 
The spectrum dN/dEQnE) is symmetric about ln£o, because 
Emax/Eo = £ , ()/£ m j n . The observed spectrum from a source 
will have this property if the pion distribution is isotropic, true 
for both WIMP annihilation and cosmic ray interactions. 

3 There is a range of predictions for the annihilation rate at the GC. 
Baryons are crucial: a massive bar can disrupt the central cusp, adiabatic 
contraction can strengthen it. Satellites are simpler as their baryons can not 
cool. The GC brightness is thus decoupled from the brightness of satellites. 
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The photon spectrum from a single pion energy must be 
convolved with the pion spectrum. We consider the process 
XX ~~ > bb, the annihilations of pairs of self-conjugate dark 
matter particles to pairs of b quarks. The process is non- 
relativistic, giving monochromatic quarks with energy E q = 
m x . The quarks each form "jets" dominated by ir mesons. In 
figure [T] we plot the phot on spectra from ann ihilations as cal- 
culated by DarkS USY dGondolo et alj|20 04). which uses re- 
sults from Pythia dSiostrand et al. 2006). The spectrum is uni- 
versal: even XX ~ * W + W~ or Z°Z° gives simila r results. Only 
the ~y~X — » r + r~ channel d i ffers appreciably dCesarini et al.l 
120041: iFomengo et alj|2004t iHooper & Tavloril2006l) . but this 
is difficult to arrange for WIMPs from supersymmetry. 

In figure [T] we show the photon spectra from sev- 
eral power- law proton sourc es as calculated by GALPROP 
dStrong & Moskalenkolll998l) . This illustrates that the power- 
law proton beams typical of astrophysical sources can not eas- 
ily mimic the gamma ray spectrum from WIMP annihilations. 
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FIG. 1 . — Spectrum of photons from hadronic processes. Solid lines depict 
annihilations to b quark pairs for several WIMP masses. The peak in the 
spectrum occurs at an energy of Epeak Kb m x /25. At low photon energies the 
spectrum is nearly independent of WIMP mass in both shape and magnitude. 
Dotted lines depict pp interactions for several proton spectra. The pp photon 
spectra are normalized to be equal at an energy of Eq = 67.5 MeV (vertical 
dotted line), where their spectral slopes are guaranteed to be equal. 



low 10 Mq due to computational limitations. The spatial dis- 
tribution of satellites is roughly spherically symmetric about 
the galactic center and extends well beyond the solar orbit, 
thus the dark matter satellites are located mostly out of the 
galactic plane. Indiv idual sources have NFW density profiles 
dNavarro et al.lll997l) , with central r~ l c usps. Satellites with 
steeper profiles, e.g. dMoore etai1ll999l) . would be easier to 
detect. We find that the brightest sources have masses in the 
10 6 - 10 7 M Q range. These brightest sources have tidal radii 
of order 100 pc, typically corresponding to 1° on the sky. We 
note that most of these objects are severely stripped. They 
have scale radii r s much larger than their tidal truncation radii 
r,, thus they have nearly pure r~ l density profiles out to r = r,. 

The surface brightne ss in gamma rays is pro portional to the 
parameter / oc J p 2 dr dBergstrom et all 1998). For a stripped 
NFW clump, at a fixed angular distance from its center, J oc 
M 2 /rf /D, where D is the distance. If the mass spectrum of 
clumps is dN/dlnM oc M~ a and the tidal radius r, oc 
the surface brightness of the nearest clump (D cx M"' 3 ) is 
/ oc M 2 " 4 ^""/ 3 . Our simulations indicate that awl and 
/3 w 1/2, thus J oc M -1 / 3 . Lower mass clumps are brighter, 
but they appear smaller (9 oc r,/D oc M^~ a ^ oc M 1 / 6 ). Less 
massive clumps would be seen as point sources. These results 
are sensitive to a and f3, thus extrapolations are difficult. 
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3. DETECTABILITY OF GALACTIC DARK MATTER SATELLITES 

The GLAST satellite d AtwoodJ [19941: iBlooml [l996t 
iGehrel s & Mic helsonl 119991) is well suited to measuring 
gamma rays from dark matter annihilations. It has an effec- 
tive area of wl m 2 , a solid angle acceptance of m 3 sr, and a 
point spread function (PSF) of 0.4° at 1 GeV energy. It will 
measure gamma ray energies between 20 MeV and 300 GeV. 
Most data will be taken in survey mode, mapping the sky with 
equal coverage with a large duty cycle. The exposure to any 
source will reach roughly 3 x 10 11 cm 2 s in a 5 year mission. 

We have estimated the number of Milky Way dark mat- 
ter satellites observable by GLAST. The calc ulation was per- 
formed with the semi-analytic method of iTavlor & Babull 
(2004, 2005a,b). The satellite mass distribution has the ex- 
pected dN/dM oc M~ 2 dGhigna et al.lll998l) . cutting off be- 



FlG. 2. — The number of detectable clumps for several WIMP masses is 
plotted against the detection threshold for counts above 1 GeV. Distinguishing 
the WIMP spectrum at 3a from a molecular cloud and a pulsar is possible at 
detection significances of 15<r (290 counts, 375 background) and 22cr (425 
counts) respectively (vertical dotted lines). 



As a fiducial case, we assume a WIMP mass of 100 
GeV and an annihilation cross-section to bb of (rxv) = 1 .6 x 
10~ 26 cm 3 s" 1 , giving 14.2 photons per annihilation above 1 
GeV Assuming a 5 year GLAST mission, and integrating 
a 1° radius around the source, the number of background 
counts is 375 (based on the EGRET extragalactic background 
dSreekumar et al.l fl998l) ). The typical brightest clump has 
(/)(1° radius) = 1400. For an example of such an object take 
2 x 1O 6 M , 3 kpc distant, tidal radius 50 pc, thus subtend- 
ing 1 on the sky. The number of signal counts above 1 GeV 
within 1° for 30, 60, 100, 150, 200, 300 GeV WIMPs is 3450, 
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FIG. 3. — Annihilation of 30 GeV WIMPs. The solid line represents the 
WIMP spectrum from a single clump with (J) = 1400 within 1°, including 
the extragalactic background. Shaded boxes illustrate the la (poisson) flux 
errors for 0.2 decade bins (the actual resolution is better than 10%, and we ex- 
pect very little dispersion between bins this large). The best fit proton power 
law (long dashes) and slope -4/3 pulsar spectra (short dashes) are indicated. 
These fits are unacceptable, including a pure power law (not shown). This 
source would have had ~50 counts above 100 MeV in th e third EGRET cat- 
alog, just above the detection limit (Hartman et al. 1999). 

1680, 900, 520, 345, 190 respectively. The number of de- 
tectable dark matter satellites is shown in figure [2] 

4. ASTROPHYSICAL SOURCES 

A pure dark matter galactic satellite has four distinguish- 
ing characteristics: hadronic spectrum from monochromatic 
quarks, spatial extent, lack of variability, and no emission at 
other wavelength s except for very diffuse inverse Compton 
and synchrotron (Baltz & Wai 2004). We will focus on the 
energy spectrum, but we note that a satellite with an NFW 
profile has a brightness proportional to 1/r, meaning equal 
flux in equal width annuli. With the 0.4° PSF of GLAST 
above 1 GeV, and the typical 1 size, we expect that the spa- 
tial extent should be detectable. Naively, emission between 
r = 0.5° and r = 1° has 58% of the total significance. 

To mimic a galactic satellite, a source would need to have 
a broken power law spectrum, no counterparts in other wave- 
lengths, and degree-scale extended emission constant in time. 
Each of these is exhibited by known sources, but none exhibits 
them all. A class of such sources would need to have identical 
spectra (as measured), as dark matter satellites would. 

In figures [3]|5] we plot the spectrum of the typical brightest 
clump ((/) = 1400) with 30, 100, and 200 GeV WIMPs to- 
gether with fits for several astrophysical source classes. The 
la errors for GLAST are shown as shaded boxes. In figure|2] 
we illustrate the detection significance required to distinguish 
dark matter from molecular clouds and pulsars. 

The gamma ray spectrum from molecular clouds is gener- 
ated by cosmic ray protons, which have a featureless power- 
law spectrum over many decades in energy. This is are exactly 
what is plotted with dotted lines in figure Q] The long dashed 
lines in figures [3]|5] show the best fit molecular cloud spec- 
tra, ruled out at high confidence in each case. Gamma rays 
from molecular clouds are expected to be extended and non- 
variable; it would be comforting to rule out counterparts, e.g. 



FIG. 4. — Annihilation of 100 GeV WIMPs. The curves are the same as in 
figure|3] Again, none of the fits are acceptable, including a pure power law. 
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FIG. 5. — Annihilation of 200 GeV WIMPs. The curves are the same as in 
figurefS] In this case, the pulsar fit is allowed at the 8% level (by simple x 2 ). 
but the others are unacceptable, including a pure power law. 

CO emission. In fact, GLAS T is likely to detec t extended, 
high-latitude molecular clouds dTorres et al. 2005]). 

Gamma ray pulsars are potentially the most problematic of 
the astrophysical sources. Their s pectra can be paramet erized 
as dN/dE cxE- r exp[-(E/E c ) a ] dNel & de Jagerll 19951) . The 
few known examples have V > 4/3. In fact, most models 
for gamma ray pulsars require this (e.g. the outer gap model 
(Romani 1996)), but V -> 2/3 is in principle possible. 

The short dashed lines in figures [3]|5] show the best fit 
gamma ray pulsar spectra. The (T = 4/3) spectrum is ruled 
out except in the 200 GeV case, where it is consistent at the 
8% level. If the low energy slope r — ► 1, the spectra become 
nearly impossible to disentangle for any WIMP mass. 

Gamma ray pulsars tend to have multi-wavelength counter- 
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parts and also tend to be near the galactic plane. A notable ex- 
ception is 3EG J 1835+59 18 which is lo cated at high latitude , 
but does have a faint X-ray counterpart (Hal pern. et.al.ll2002h . 
The well known radio quiet gamma ray pulsar Geminga is lo- 
cated within 5° of the galactic plane. 

The variability of the pulsar is difficult to determine in a 
blind search of the period-period derivative plane. To mimic 
an extended galactic satellite, a cluster of pulsars would be 
required, with no counterparts in other wavelengths. 

Plerions will typically have multi-wavelength counterparts, 
especially in X-rays, and are located close to the galactic 
plane. They are compact sources in X-rays (~ 1'), but at 
GLAST energies they may be detected as extended sources. 

Supernova remnants will have power-law spectra, convinc- 
ingly ruled out. They will also have multi-wavelength coun- 
terparts and are likely to be near the galactic plane. 

Blazars are variable point sources with a power-law gamma 
ray spectrum, and have counterparts. They contradict all four 
necessary qualities of dark matter satellite emission. 

5. DISCUSSION 

We have shown that the brightest dark matter satellites seen 
by GLAST should be distinguishable from other known types 
of astrophysical sources, for WIMPs less massive than about 
150 GeV and cross sections of (ctv) = 1.6 x 10~ 26 cm 3 s" 1 . 
Sub-s ubstructure could en hance these signals by a factor of 
a few (Diemand et al. 2007), extending the accessible masses 
and cross sections. Dimmer sources (extended only, to min- 
imize pulsar contamination) could be stacked, improving the 
discrimination. Any sources mimicking dark matter annihi- 
lation would be very interesting, and would be compelling 
targets for study in a multiwavelength campaign. 

Dark matter sources are excellent targets for Imaging At- 



mospheric Cerenkov Telescopes (IACTs). The mass of the 
WIMP must be above the IACT analysis threshold, at the 
present time around 100 GeV. The sensitivity of IACTs is cur- 
rently limited by the charged particle background. 

A follow-up campaign of 500 hours with an IACT of 0.2 
km 2 on our brightest clump, taking the 100 GeV model, 
can provide a 5rr detection of line emission from the pro- 
cess XX ~ y 77' f° r a branching ratio of B = 1.2%. This as- 
sumes 99% rejection of hadrons and 15% energy resolution. 
If the hadron rejection were improved by a factor of 10, the 
electron background dominates and line sensitivity improves 
to B = 0.005. If the electron background were also elimi- 
nated, the extragalactic gamma-ray background would limit 
the sensitivity to B = 0.0003. The predicted branching ratio is 
B ~ 0.001. Obviously, the detection of a line at these energies 
would demonstrate the existence of particle dark matter. 

The Large Hadron Collider (LHC) may discover a can- 
didate WIMP, and measure its mass at the 10% level on 
a timescale that matches the GLAST program. A simple 
estimate shows that GLAST can constrain the mass at the 
25% level, for a 100 GeV WIMP. If the GLAST and LHC 
mass estimates match, the WIMP hypothesis would be greatly 
strengthened. With strong evidence for particle dark mat- 
ter in hand, including accelerator measurements of cross sec- 
tions dBaltz et al.ll2006l) . it would become possible to map the 
galactic dark matter in the gamma ray sky. 
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